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We consider the two-dimensional Lorentz gas with Poisson distributed hard disk scatterers and 
a constant magnetic field perpendicular to the plane of motion. The velocity autocorrelation is 
computed numerically over the full range of densities and magnetic fields with particular attention 
to the percolation threshold between hopping transport and pure edge currents. The Ohmic and 
Hall conductance are compared with mode-coupling theory and a recent generalized kinetic equation 



valid for low densities and small fields, 
magnetic field. 



We argue that the long time tail as t persists for non-zero 
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I. INTRODUCTION 

Two-dimensional electron films can be manufactured 
with high perfection in GaAs heterostructures. At low 
temperatures a mean free path of over 10 4 nm is reached 
and to a very good approximation the electrons may be 
considered as noninteracting. To have some interesting 
physics one nanostructures the probe by lithographic or 
other techniques. Thereby a strongly repulsive potential 
is imposed on the electrons with a maximum above the 
Fermi energy (quantum antidots). If the imprinted struc- 
ture is on a scale larger than the Fermi wave length, ad- 
justable to be of the order of 50nm, one hopes to capture 
the transport properties already in the classical approxi- 
mation. 

So far, the most popular geometry has been a regu- 
lar array of antidots which corresponds classically to the 
Sinai billard and should result quantum mechanically in 
the Hofstadter butterfly. The magneto-transport of this 
periodic structure has been studied in great detail, both 
experimentally and theoretically In our note 

we investigate randomly placed antidots. To our knowl- 
edge, the so far best experimental realization has been 
achieved by Lutjering |)|. We compare our results with 
his measurements in the conclusions. 

In kinetic theory randomly distributed scatterers are 
known as the Lorentz gas, which has proved to be an 
important testing ground. In particular, one can under- 
stand precisely the assumptions for the validity of the 
(linear) Boltzmann equation |l0[ and check on the accu- 
racy of the low density expansion and its non-analytic 
character []llj-|l6f . Also the long time tails in the velocity 
autocorrelation function are seen most convincingly in 
the Lorentz gas (l7]]. In distinction from the work men- 
tioned we investigate here the dynamics in presence of a 
magnetic field perpendicular to the plane of motion. 

The model has a strong geometric flavour: one places 
randomly disks (the scatterers) of radius a in the plane 
at density n s . The disks may overlap. In the region out- 
side the disks we have independent point particles with 



density n e . They have mass to*, charge e, and move in a 
uniform external magnetic field B ex . Thus a single parti- 
cle travels along a circle and is elastically reflected upon 
collision with a scatterer. We denote the velocity of the 
particle at time t by u(i). Clearly |u(t)| is conserved and 
we set it equal to the Fermi velocity vp , since at low 
temperatures contributions to the transport only come 
from the Fermi surface. The radius of gyration is then 
R e = v F m*/eB ex . 

We are interested in the magneto-transport which re- 
lates the steady state current j to an in plane uniform 
electric field E by 



j = aE 



(1.1) 



for small E and want to understand how a depends on 
B ex and n s . a±i and 022 are the Ohmic conductivities. 
In our case an = 022 by isotropy. a 12 — —&21 is the Hall 
conductivity. The magneto-transport will be studied in 
linear response. This means the dynamics is the one just 
explained (zero electric field) and the transport coeffi- 
cients are given in terms of the time-integrated velocity 
autocorrelation functions. Physically one has to average 
over all of phase space. For a large sample, this is equiv- 
alent to fixing the initial position and averaging over the 
scatterer distribution. It is this prescription we will use 
both theoretically and in the numerics. Before spelling 
out the details we introduce dimensionless quantities. 

Space is measured in units of the disk radius a and 
velocity in units of the Fermi velocity vp. The scat- 
terers have then radius 1 and their dimensionless den- 



sity is p 



The dimensionless radius of gyration 



is R = R e /a and the corresponding magnetic field is 
B = 1/R = eaB ex /vFm* . Let v(t) — u(t)/v F be the 
velocity of the particle at time t starting at the origin. 
Clearly \v(t)\ — 1. v(t) depends on v(0) and the partic- 
ular configuration of the scatterers. 

We define the dimensionless velocity autocorrelation 
function by 



Cij(t) = (Vi(p) Vj (t)) ,i,j = 1,2. 



(1.2) 



Here ( • ) is a double average. Firstly there is an av- 
erage over scatterers. The centers of the disks are dis- 
tributed according to a Poisson process with uniform 
density p conditioned on the set {x | |x| < 1} being free 
of centers. Secondly we average over the initial velocity 
v(0) = (cos tp, sin ip) uniformly in (p. (By rotational in- 
variance this second average could be omitted, but it is of 
advantage numerically.) The conductivity tensor is then 



n e e 



*(ui(0K(i)). (1.3) 



D depends on p and B. Also of interest is the frequency 
dependent conductivity defined by 



(Tij{u) = -V AiH , D 13 (uj) = / dte wt { Vl {Q) Vj {t)) 



(1.4) 

We note that (vi(0)vj(t)) — (vi(0)vj(—t)) by stationar- 
ity and («i(0)ua(t)) = -(«i(0)«2(-t)) = -(w 2 (0)«i(t)) 
by time reversal. Therefore D\2 = —D%x and Da = 

(l/4)/-"<ft(v(0)-v(t)>. 

As stated in ( ]1 .3| ) the conductivity is ill-defined. There 
is always a non-zero probability that the particle will not 
be scattered at all. If so, (vi(0)vi(i)) = (1/2) cos(Bi) 
and (vi(0)v2(t)) = (1/2) sin(Bi), and the time integral 
in (1.3) needs a reinterpretation. Physically there will al- 
ways be a weak elastic scattering by impurities, i.e. once 
in a while the velocity direction is randomized. In ap- 
proximation the velocity autocorrelation function is then 
modified to 



Thus for the Lorentz gas there should be a well-defined 
steady state current j for B ^ 0. We did not find a gen- 
eral argument to establish its existence. The situation at 
B = is very different. Then the particle is accelerated 
along the direction of E. Since only the velocity direc- 
tion is randomized by collisions, the energy input from 
the electric field cannot be dissipated and no meaningful 
steady state current is reached as t — > oo. However for 
small E the time dependent current settles at a plateau 
over a time span the lon ger the smaller E, whose value 
equals ctE with a of ( |l.3| ). 

To give a brief outline: In the following section we dis- 
cuss the dependence of D on p, B with particular atten- 
tion to the two percolation thresholds. We also compare 
our numerics with the mode coupling theory of Gotze and 
Leutheufier which seems to be the only theoretical 
prediction at intermediate densities. Recently Bobylev et 
al. |p^|| derived a generalized transport equation for low 
scatterer densities and small fields. In the appropriate 
domain of validity their predictions are in fact very accu- 
rate (Section 3), and clearly improve on the phenomeno- 
logical Boltzmann equation with magnetic field. One of 
the famous results on the Lorentz gas at B = is the slow 
decay of the velocity autocorrelation function as — t~ 2 for 
large t in two dimensions |2(J. For densities p < 0.25 
such a power law has been well established numerically 
p7| . At larger densities there is a pre-asymptotic decay 
approximately as —t~ 1A and, with reasonable numerical 
effort, the true asymptotics cannot be seen any more. In 
Section 4 we argue that also for B / the velocity au- 
tocorrelation decays as t~ 2 for large t. We conclude by a 
comparison with the experiments of Lutjering and with 
some comments. 
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Vi(0)v j(t)) 



and the proper definition reads 
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Dij(u>) = lim / dt 
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Vi(0) Vj (t)) 



At zero density, no scattering, (|1.6|) results in 
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(1.5) 



(1.6) 



as well known from the Drude theory of the Hall effect. 
In particular, at lu = both the Ohmic conductivity and 
the Ohmic resistance vanish. 

It is instructive to go back for a moment to the case 
of an in plane electric field E. One finds that for freely 
moving particles the average current approaches 



— ( — Ei, Ei ) 

2B \ 2, LJ 



(1.8) 



as t — * oo in accordance with (1.7) at u> = 0. In fact, one 
would expect that scattering cannot make things worse. 



II. THE STATIC CONDUCTIVITY 

To discuss D in its dependence on p, B at zero fre- 
quency it is useful to consider first the limiting cases. As 
explained D 11 (0,B) = 0, D 12 (0,B) = 1/2B at p = 0. 
On the other hand for B = we have D^p, 0) = 0. 
D\\(p, 0) has been studied numerically ■ For p — > 

one obtains the Boltzmann value Dn(p) = 3/16/5. Thus 
close to (p, B) = 0, D is somewhat singular and roughly 
of the form (p 2 + B 2 )^ 1 / 2 . Its precise functional depen- 
dence will be discussed in Section 3. 

As the density is increased the disks percolate. This 
means that for p > p c , p c = 0.36 with probability one 
the origin is contained in a finite domain bounded by 
scatterers. For p < p c the origin is connected to infinity 
by a path not intersecting the scatterers. Close to p c 



11.5 
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numerically one finds Dn(p, 0) = \p — p, 
with no theoretical explanation yet. 

In fact, the Lorentz gas has a second percolation 
threshold. Let us fix p < p c and increase B. Above 
some critical value B c , the trajectory will either be a 
circle or skip along a , possibly large, cluster of a finite 



number of disks. The bulk current for B < B c is reduced 
to a pure edge current. Since the particle cannot leave a 
cluster, the mean square displacement is bounded which 
implies Dn(p,B) — for B > B c by the Einstein rela- 
tion 2Dn = lim 4 ^ 00 (r 2 (t)) /2t. A hopping of the particle 
from cluster to cluster is only possible, if disks with radius 
1+R percolate which means, in our units, p(l + i?) 2 > p c , 



Vpc/p- i 



(2.1) 



Strictly speaking the B c of ( |2.1|) is only an upper bound 
on the true B c . One could imagine that already for a 
slightly smaller B hopping is suppressed. Numerically, 
we see a smooth variation through B c and such a fine 
point cannot be decided. In Fig. 1 we plot the two 
domains in which the Ohmic conductivity vanishes. 

The behavior of Dyi is somewhat more complicated. 
We always have the contribution of the circle orbits. Ac- 
cording to the Poisson distribution their probability is 
exp[-7rpK(i?)] with k(R) = R(R + 2) for R < 1 and 
k(R) = 4R — 1 for R > 1. We decompose the average 
correspondingly as 



{vi(0)vj{t))>i + (MO)vj(t))o-. 



(2.2) 



where ( • )>i is the average over all initial v(0) and all 
scatterer configurations such that there is at least one 
collision for < t < oo. The normalization is = 1 — 
exp[— Trpn(R)]. If (vi(0)vj(t))>i is absolutely integrable, 
then 

Da = lim (vi(0)[ Xj (t) - or i (0)])> 1 + (1 - S^-e-*^™ , 

t — >oc Z 

(2.3) 

where x(£) = L ds v(s) is the position of the particle at 
time t. The average (vi(0)[xj(t) — Xj(Q)])>i can be fur- 
ther decomposed into averages counting with how many 
different disks the particle collides in the course of time. 
The lowest contribution corresponds to the particle skip- 
ping around a single disk, etc.. Now, if either B > B c 
or p > p c , then by assumption the sum over all clus- 
ters with a finite number of disks exhausts already the 
complete average (vi(0)[xj(t) — Xj(0)])>i for arbitrary t. 
(In the domain where Dn > 0, the particle collides with 
an infinite number of disks). Let us fix then a particu- 
lar finite cluster of disks and let T be the set {x, v} of 
initial conditions such that the particle will eventually 
collide with each one of the disks in the cluster and no 
others. By assumption |x| < const, for {x, v} G T. If 
the dynamics in T is mixing (or decomposes into mix- 
ing components), then Um t _ KXJ («x(0)[a;2(i) — X2(0)])r = 
<«i(0))r<sa(0))r-(«i(0)x 2 (0))r = -(«i(0)x2(0))r, since 
(«i(0))r = (-3j x i{t))r = by stationarity. 

For one disk the motion is integrable [glj and the mix- 
ing assumption fails. In the Appendix we compute the 
contribution D± 2 , to ( |2.3| ) coming from rosette orbits 



around a single disk. The case of two non-overlapping 
disks was studied in J21J] . They still found elliptic islands, 
however with a small measure already. The mixing as- 
sumption seems to bc satisfied approximately. One might 
hope that the remainig contribution, — (i>i(0)x2(0))r, in- 
tegrated over allowed scatterer configurations averages to 
zero. Numerically this does not seem to be the case. At 
density p = 0.1, the circling and rosette orbits together 
account only for 50% to the observed -D12 in the range 
1.5 < B < 2.5, cf. Fig. 2b. 

To have a more complete picture of D we simulate 
the Lorentz gas numerically. For given B and scatterer 
configuration we compute Vj(t) up to 60 collision times. 

The system size is chosen so large that x(i) = f Q ds v(s) 
never hits the boundary. To speed up the simulation we 
use a hierarchical search for the next point of collision. 
For each scatterer configuration we average over 100 ran- 
domly chosen initial velocity directions. For Cy(i) to 
be sufficiently smooth typically one has to average then 
over 1 6 s ample paths. The conductivity is determined 
from ( [2.3[ ). In most of parameter space (vi(0)xj(t)} has 
not yet reached its asymptotic value, which reflects the 
slow decay of the velocity autocorrelation functions. We 
essentially extrapolate "by hand" to t — + 00 which re- 
sults in a slight overestimate whenever there is an in- 
dependent check. In Figs. 2-5 we display our results 
at densities p = 0.1, 0.15, 0.2 and 0.3. The percolation 
threshold is indicated by a vertical line. Note that for 
p = 0.1 the B-scale starts at B = 0.5. Our data show a 
fairly smooth interpolation of the asymptotics at B = 0, 
B = 00. The most surprising feature is an initial increase 
of the Ohmic conductivity with B at intermediate den- 
sities. Apparently, the curved trajectory can bend itself 
more easily through the dense "labyrinth" of scatterers. 
The Hall conductuctivity rises steeply to its maximum 
and then levels off. The maximum is shifted to smaller 
B as the density decreases. In Fig. 2b we ins ert t he con- 
tribution of the circling orbits according to (2.3) and of 
rosette orbits according to (AE). 

To our knowledge the only attempt to derive the 
magneto-conductivity at intermediate densities is the 
mode coupling theory of Gotze and LeutheuBer pl|. In 
Fig. 6 we compare their prediction with our simulation 
data, note the particular choice of units. For the low- 
est density, p — 0.032, we find good agreement. Essen- 
tially the same behaviour is obtained from the general- 
ized Boltzmann equation, to be discussed in the following 
section. However, the kinetic theory fails already at den- 
sity p — 0.13 whereas mode coupling is still a reasonable 
approximation. The next higher density p — 0.25 can no 
longer be accounted for. 



III. LOW DENSITY, SMALL FIELDS 

For B = the linear Boltzmann equation is exact at 
low density. More precisely, for p — > collisions of the 



particle with the same scatterer become unlikely and the 
density of particles in phase space on the scale of the 
mean free path is governed by the Boltzmann equation. 
In particular this yields Du(p,0) = 3/16/9 as p — > 0. 
With external magnetic field B one generalizes in the 
obvious way to 



<9 t /(x, <f, t) 
+2p 



= (— cos tp di — sin ip d 2 — B d v )f(x., (p, t) 

V I sin |y|[/(x, <p - <p', t) - /(x, <p, t)] . (3.1) 



Here / is the distribution function at x, 

v = (cos ip, sin ip) , t. On the basis of (3.1) one obtains 



-^11 



1 



8p/3 



B 



2 (Sp/3) 2 + B 2 ' LZ 2 (8p/3) 2 + B 2 
which we rewrite in scaling form as 

1 



(3.2) 



u \\ — 



B 2 



1 



B 2 



gUp/B), 



9°u(p/B), 



8y/3 



(8y/3) 2 
1 



(3.3) 



For a rigorous derivation the radius of gyration must 
be of the order of the mean free path. Thus B must 
vanish linearly with p. As observed by Bobylev et al. 
fl9f even in the limit p — > some recollisions survive. 
This is most easily seen for circle orbits with no col lisions 
at all. According to the Boltzmann equation (3A) even 
after several turns the particle still has some probability 
to be scattered. However, for the mechanical Lorentz 
gas after one completed revolution the annulus is surely 
free of disks and no scattering events can occur. In ]19| ] 
the circling orbits and the recollisions are properly taken 
into account for p — > and a generalized kinetic equation 
with memory term is derived. On this basis the velocity 
autocorrelations arc com puted. The conductivity is still 
of the scaling form (fo.Sl), 



D* — 



1 



D* 



y/p 2 + B 2 

1 

yjp 2 + B 2 



9u(p/B) 



9*MB) 



(3.4) 



which is just a consequence of B/p = const, as p - 
However the scaling functions are now modified to 



9n(y) 



1 



y 



(2|ry(a0) 5 



-(l-x 2 )2yj(x), 



(3.5) 
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(l + x 2 (2 yl (x)) 2 ) 



x = c- 2 ^ . 



j(x) = 1 - 
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(3.6) 



(3.7) 



(3.8) 



In Fig. 7 we plot , g[ 2 . Note that in terms of the polar 
angle in the p-jB-plane the scale is highly compressed at 
the left. In Fig. 8 we plot the correction to the Boltz- 
mann value ( |3.3| ). 

Of course, D* reproduces the correct limiting be- 
haviour for p — > and B — * 0. According to ( pTl| ) the 
percolation boundary behaves as B c = l.QQ^fp for small 
p. Since for the validity of (3.4) B is scaled proportional 
to p, D* cannot see this threshold. 

In Figs. 9, 10 we compare our numerical results with 
D*. We also include D°. As expected, D* is a con- 
siderable improvement. Note that the agreement is not 
uniform in y 1 B 2 + p 2 . This reflects that at B = the 
Boltzmann equation has a restricted range of validity, 
e.g. D u (p,0)/ J D? 1 (p) = 0-65atp = 0.1. 



IV. VELOCITY AUTOCORRELATIONS 



On the level of the linear Boltzmann equation (3.1) one 
obtains an exponential decay for the velocity autocorre- 
lations. Explicitely 

C°n(t) = ie-l*l/^ cos(Bt) , C° 12 (t) = Vl'^ sm(Bt) , 

(4.1) 

with To = (8p/3) _1 . The correct low density/ small field 
behaviour has a more interesting structure. To state the 
result it is convenient to introduce the Laplace transforms 



F(z)= dte- zt 2(Cn(t)+iC 12 (t)) 



(4.2) 



for ^te{z) > 0. In the same approximation as leading to 
.41) one obtains 



F(z) 



1 _ c^+^r 



-(z+v)T 



z-iu + v j #ff('0) 1 _ e( :L z+ !,')T+i 



(4.3) 

We find it convenient to compare the Fourier transforms 



(4.4) 



with the simulation data. Since Cu(t) is even and C\2{t) 
is odd, we have Cn(ui) = $te(F(iui) + F(—iuj)), Ci2(o>) = 
idie(—F(iuj) + F(—iuj)). The circle orbits yield 5-peaks 
at lo — ±1/R with weigh t ex p(— 4wpR). In Fig. 11 we 
plot the prediction from ( |4.3| ) and compare it with the 
numerics. Note that the 5-peaks are out of scale. As 
anticipated the agreement is excellent, in fact over the 
whole low density/ small field regime. We plot 612(0;), 
which turns out to be negative, Fig. 12. As the ratio 
p/B is increased, the characteristic double peak merges 
into a single peak which shifts then to u> = 0. 

Away from low density/ small fields we have no theory 
to compare with. For B = Cu(t) has the long time tail 
of the form C n (t) = -at- 2 with a > 0. We refer to |J] 
for a detailed discussion. The heuristic explanation runs 
as follows: At long times the main contribution of Cn(t) 
comes from paths returning to the origin at time t. Let 
v(0) = (1,0). Then the particle will be more likely to 
return from the right which yields a negative correlation 
in the velocity. For the excursion away from the origin 
we use a random walk approximation. If at some inter- 
mediate time the particle arrives at the line x = 0, then 
it will return to the origin equally likely from right and 
left and the contribution to Cn(t) vanishes. Therefore 
the tail can be computed from a return to the origin at 
time t of a random walker without ever hitting the line 
x = 0. This probability decays as t~ 2 in two dimensions. 
Clearly, in our argument we only used that the motion 
before returning to the origin is diffusive. This remains 
valid for non-zero B. Thus Cu(t) and C\2(t) should have 
a decay as t~ 2 for t — > 00. For B = the long time tail 
is most clearly seen at p = 0.15 fl7|| . We increase B to 
0.2 and average over 3 • 10 7 sample paths. The effective 
exponent for both correlations is 1.9 with a negative pref- 
actor, Fig. 13. However the prefactor of Cn(t) becomes 
positive in the range B — 0.5 ... 0.8 and of C\2{t) in the 
range B = 0.6. ..3. Thus when the particle returns to 
the origin it picks up more complicated velocity correla- 
tions than in the case B = 0. As one example of such 
a sign reversal in the 1, 2-correlation we display the data 
for p = 0.2, B = 2.3, Fig. 14. The effective exponents 
are approximately 1.4. 

V. DISCUSSION 

We studied transport in the two-dimensional Lorentz 
gas with a constant magnetic field perpendicular to the 
plane of motion. The theory of Bobylev et al. is in fair 
agreement with our simulation data, both for the trans- 
port coefficients and the velocity autocorrelations. Away 
from low density/small magnetic fields there is little the- 
ory to compare with. The qualitative properties of the 
magnetotransport can be guessed from the well under- 
stood limiting cases B — > and B — > 00. The most 
unexpected feature is an increase in the Ohmic conduc- 
tivity with increasing B at intermediate densities. Such 



a behavior has also been found experimentally Q. 

In these measurements antidots are imprinted at ran- 
dom locations with densities l/(1000nm) 2 , l/(600nm) 2 , 
l/(400nm) 2 , l/(300nm) 2 , and l/(240mn) 2 , resp. . 
The electrons feel a screened potential. For the pe- 
riodic case this is usually modelled by V(xi,X2) — 
Vo\sm(irxi/a)sm(TTX2/a)\P with (3 ranging from 2 to 4. 
For a random distribution the true potential is more com- 
plicated, in particular when there is strong overlap. Thus 
we cannot expect quantitative agreement between our 
hard disk model potential and the experiment Q . There 
the reduced density roughly ranges from 0.01 to 0.25 in 
our units. The dimensionless magnetic field varies from 
to 3. Qualitatively the conductivities D\\, D12 in de- 
pendence on B follow our curves. However the measured 
D\2 is smaller by approximately a factor 1/2. At the 
largest density a distinct increase in D\\ with B close 
to B = is observed. As in our numerical studies the 
percolation threshold, B c , is hardly seen in the Hall con- 
ductance. The Ohmic conductivity is small for B > B c , 
but the transition region is fairly broad. 

The velocity autocorrelations have a slow decay over 
the full range of parameters. Presumably it is governed 
by i -2 , which is however severely masked by an even 
slower preasymptotic decay. The definite sign of the pref- 
actor at B — is not retained. 
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APPENDIX A: HALL CONDUCTIVITY OF 
ROSETTE ORBITS 

We determine the contribution to D12 from orbits 
which circle around a single scatterer, cf. Section 2. Since 
B > B c , we discuss only the case R < 1. There is a cor- 
responding formula for R > 1. 

We choose a large square box A. Let y — {yi, ...,y„} 
be the centers of the scatterers and P p (dy) their nor- 
malized Poisson distribution at density p. We set 
A(y) = {x e A I |x-y 3 -| > 1, j = l,...,n} and de- 
note by |A(y)| the area of this set. Let Fj(t;x,ip,y) — 
cos if sin ip(t; x, <p, y) if ip(t) defines a rosette orbit at scat- 
terer j and Fj = otherwise. Then 

(«i(0)«a(t))i = J Ppi^J]^ J J ^E^W- 

My) ' j=1 

(Al) 

We have |A(y)| = |A|e~ 7rp . By translation symmetry we 
shift the scatterer j to the origin and introduce the po- 
lar coordinates x = (r cos a, r sin a). Let f r (t;a,tp) — 



cos ip sin (f{t; r, a, tp) if the orbit with initial center of gy- 
ration (r cos ev, r sin a) and initial velocity (cos ip, sin ip) 
defines a rosette orbit around the scatterer located at 
and f r — otherwise. Let d{y) — min,,- | y ^ | . Taking the 
limit |A| — > oo we obtain 

(vi(0)v 2 (t)) 1 = pc^J P p (dyM2 < d(y)} x 

(l+R)A(d(y)-l-R) 

drr— da d(pf r (t,a,ip) = 



l-R 



pe*P J P p {dy) X {r + R+l< d(y)} x 

l-R 



l+R 



pe" 



* p f drre-" (r+it+1)2p -^ f da f ckpf r (t,a,<p) (A2) 



2tt 



l-R 



with x(A) denoting the characteristic function of the set 
{A}. 

As before, the Hall conductivity of rosette orbits is 
given by 



D$ = lim / dte~ zt (v 1 (0)v 2 (t)) 1 



(A3) 



By rotation invariance, the a integration becomes trivial. 
Thus we have to determine 



lim / dte~ zt 

z^O 





J d(pf r (t;a,cp) 



(A4) 



We fix a = \ and r such, that 1 - R < r < l + R. Then 



-f+ < P < V+j with 

ip+ = arccos((l - R 2 - r 2 )/(2Rr)) . 



(A5) 



(p(t) increases linearly in t. At a collision tp(t) jumps by 
the angle 



Thus 



* = 2arccos((r 2 - R 2 - 1)/2R) 



tp(t) = <p+-+n(t,<p)V, 



(A6) 



(A7) 



where n(t, tp) is the number of col lisio ns up to time t for 
the initial ip(Q) = (p. Inserting in (A4) yields 



oo V?+ 

lim J dtc~ zt J dip cosy?sin(</j + + n(t, (p) 1 ^) 



lim R / dip cos(p(2i(— zR + i)) 1 x 



-ip + 



oc 

e i{n+l)^ e (2n+l)(-zR+i) V+ e zR. V ^ +cc 



■ ■ 

i? / (it^ cos <^[(cos p — cos <^-i-) 



-<p+ 

(2(l-cos(2^+ + *)))- 1 x 

( cos(3p+ + *) - cos(< < 5 + + *))] = 

i?(<^+ - isin2^ + ) + 

2i?(sinv3 + ) 2 sin(2^ + + *)(1 - cos(2<^ + + . 

(A8) 

Combined with (|A2| ) we obtain 

1+-R 



D$=R J drrpe^ 1 - (r+R+1 ^W + - 



sin ip+ cos ip+ 



l-R 



2 (sin^+) 2 cos(^+ + (*/2))/ sin(^ + + (*/2))} 



(A9) 



In Fig. 2b we plot our result (|AJj) for p = 0.1 as a 
function of B = 1/R. 



-¥>+ 



[1] A. Lorke, J. P. Kotthaus, and K. Ploog, Phys. Rev. B 

44, 3447 (1991) 
[2] K. Ensslin and R. Schuster, in III-IV Semiconductor 

Quantum Systems, ed. by K. Ploog, Institiution of Elec- 
trical Engineers, New York, IEE 1996 
[3] D. Weiss, Phys. Rev. Lett. 66, 2790 (1991) 
[4] T. Schlosser, Quantentransport in lateralen Ubergittern, 

Dissertation, Universitat Miinchen, 1995 
[5] T. Geisel, A. Zacherl, and G. Radons, Phys. Rev. Lett. 

59, 2503 (1987); Z. Phys. B 71, 117 (1988) 
[6] R. Fleischmann, T. Geisel, and A. Ketzmerick, Phys. 

Rev. Lett. 68, 1367 (1992) 
[7] M. Fliefier, G. J. O. Schmidt, and H. Spohn, Phys. Rev. 

E 53, 5690 (1996) 
[8] R. Silberbauer, Quantenmagnetotransport und Quanten- 

chaos in Antidot-Ubergittern, Universitat Regensburg, 

1994 

[9] G. Lutjering, Nanostrukturierte Antidot-Systeme: vom 
ballistischen Billard zum Metall-Isolator Ubergang, Dis- 
sertation, Universitat Stuttgart, 1996 

[10] H. Spohn, Large Scale Dynamics of Interacting Particles, 
Berlin, Springer 1991 

[11] M. H. Ernst and A. Weijland, Phys. Rev. A 34, 39 (1971) 



[12] A. Wcijland, J. Math. Phys. 15, 1942 (1974) 

[13] W. Gotze, E. Leutheufier, and S. Yip, Phys. Rev. A 23, 

2634 (1981); A 24, 100 (1981); A 25, 533 (1982) 
[14] E. Leutheufier, Phys. Rev. A 28, 1762 (1983) 
[15] A. J. Masters and T. Keyes, Phys. Rev. A 26, 2129 

(1982); A 27, 2603 (1983) 
[16] B. J. Alder and W. E. Alley, J. Stat. Phys. 19, 341 (1978) 
[17] C. P. Lowe and A. J. Masters, Physica A 195, 149 (1993) 
[18] W. Gotze and E. Leutheufier, Z. Physik B 45, 85 (1981) 
[19] A. V. Bobylev, F. A. Maao, A. Hansen, and E. H. Hauge, 

Phys Rev. Lett. 75, 197 (1995) 
[20] B. J. Alder and T. E. Wainwright, Phys. Rev. A 1, 18 

(1970) 

[21] N. Berglund, H. Kunz, J. Stat. Phys. 83, 81 (1996) 



Figure captions 



Fig. 1: The percolation thresholds for the Lorentz gas 

Fig. 2: Ohmic conductivity Dn (a) and Hall conductivity D\ 2 (b) as functions 

of B for p = 0.1 in dimensionless units. ( ) is the zero-scatterer contribution 

of (2.3) and ( ) is the zero- and one-scatterer contribution of (A9) 

Fig. 3: Ohmic conductivity Dn (a) and Hall conductivity D\ 2 (b) as functi- 
ons of B for p = 0.15 in dimensionless units 

Fig. 4: Ohmic conductivity Dn (a) and Hall conductivity D\ 2 (b) as functi- 
ons of B for p = 0.2 in dimensionless units 

Fig. 5: Ohmic conductivity Dn (a) and Hall conductivity D\ 2 (b) as functi- 
ons of B for p = 0.3 in dimensionless units 

Fig. 6a: Ohmic conductivity Dn according to the mode-coupling theory of [18] 
for various densities 

Fig. 6b: Simulation data of the Ohmic conductivity Dn for the same densi- 
ties and in the same units as Fig. 6a 

Fig. 7: The scaling functions g n (a) and g[ 2 (b) as functions of x = e~ 2vp l B 

Fig. 8: The ratios D^/D^ (a) and D* 12 /D° 12 (b) as functions of x = e~ 2vp / B 

Fig. 9: Inverse Ohmic conductivity 1/Dn (a) and inverse Hall conductivity 

1/D 12 (b) as functions of (B 2 + p 2 ) 1 ' 2 for fixed ratio p/B = 0.025. ( ) is 

the Boltzmann theory, ( ) the improved Boltzmann theory, and (<)) are the 

simulation data 

Fig. 10: Inverse Ohmic conductivity 1/Dn (a) and inverse Hall conductivity 

1/D 12 (b) as functions of (B 2 + p 2 ) 1 ^ 2 for fixed ratio p/B = 0.1. ( ) is the 

Boltzmann theory, ( ) the improved Boltzmann theory, and (<)) are the simu- 
lation data 

Fig. 11: Real part of Cn{u) with p/B = 0.025 and (B 2 + p 2 fl 2 = 0.105 from 
the improved Boltzmann theory (a) and from the simulation data (b) 

Fig. 12: Negative imaginary part of C'i 2 ((jj) with p/B = 0.025 and (B 2 + p 2 ) 1 ^ 2 = 
0.105 from the improved Boltzmann theory (a) and from the simulation data (b) 



Fig. 13: Integrated velocity autocorrelation functions (v\(0)xi(t)) (a) and {vi(0)x 2 (t)) 
(b) with p = 0.15 and B = 0.2 as functions of (r/t) - 9 

Fig. 14: Integrated velocity autocorrelation functions (v\(0)xi(t)) (a) and {vi(0)x 2 (t)) 
(b) with p = 0.2 and B = 2.2 as functions of (r/t) 0A 
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